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Abstract 

We find quasinormal modes of near extremal black branes by solving a singular bound- 
ary value problem for the Heun equation. The corresponding eigenvalues determine the 
dispersion law for the collective excitations in the dual strongly coupled M = 4 supersym- 
metric Yang-Mills theory at finite temperature. 
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1. Introduction 

The gauge theory /gravity correspondence |1],§,[§ provides useful insights into the prop- 
erties of a strongly coupled supersymmetric Yang-Mills (SYM) theory at non-zero temper- 
ature. A well-known example is the finite-temperature M = 4 SU(N) SYM theory in 4d 
which in the large N, large 't Hooft coupling limit is dual to the gravitational background 
of N parallel near extremal three- branes, with temperature related to the parameter of 
non extremality. For this theory, a number of quantities such as the free energy |4],[5], 
the Wilson loop |||7|], the shear viscosity |@|§, the -R-charge diffusion constant ||] , and 
the Chern-Simons diffusion rate [JTD| have been computed using the methods of gauge the- 
ory/gravity duality. Since supersymmetry is broken at finite temperature and thus no 
non-renormalization theorem is expected to hold, quantities computed in the strong cou- 
pling regime using a gravity dual differ from their analogues obtained at weak coupling 
via perturbation theory. The results of a gravity calculation are then regarded as a pre- 
diction for the SYM theory, assuming that the AdS/CFT correspondence is valid at finite 
temperature^. 

Dynamical properties of a thermal gauge theory are encoded in its Green's functions. 
In the context of AdS/CFT, Minkowski space Green's functions can be computed from 
gravity using the recipe given in [[TU]. Unfortunately, for a non-extremal background, 
only approximate expressions for the correlators are usually obtained. For example, the 
retarded propagator of the gauge invariant local operator F 2 (dual to the dilaton) defined 
by 

G^^-ij^e-^ m{l F Wm (1.1) 

can be explicitly computed only at zero or very high (with respect to the absolute value 
of momentum) temperature. At zero temperature, the retarded propagator ( |1.1|) has a 
branch cut singularity for \u\ > |k|, 

<?«(„, k ) = iV2( ~ 6 ^ k2)2 (In | - u 2 + k 2 | - ine(u 2 - k 2 ) sgnu,) . (1.2) 

In the high temperature limit ui/T <C 1, |k|/T <C 1, the propagator is analytic in the 
complex w-plane, 

N 2 T 2 

G R (u, k) = (i 2nTu + k 2 ) . (1.3) 

16 



1 A non-trivial check of the validity of the correspondence at finite temperature was made 
recently in M . 



1 



However, for generic values of u and k, we expect G R (u,k) to have poles corresponding 
to the spectrum of collective excitations of the SYM plasma. 

One can compare the situation to the simpler case of the 2d CFT dual to the Bahados- 
Teitelboim-Zanelli (BTZ) black hole background. There, the retarded Green's functions 
can be computed exactly. For illustration, consider the case of the conformal dimension 
A = 2. Then 



G* d (u,k) 



u 2 - k 2 

471-2 



u — k\ ( . uj + k 



4 ,T 



;i.4> 



where we have put Tl, = Tr and ignored the constant prefactor for simplicity. The high 
temperature limit of ( |1.4j) is an analytic function of u. In general, however, G^u, k) has 
infinitely many poles located at 

uj n = ±k-i4nT(n + l), n = 0, 1, . . . . (1.5) 

When T — > 0, the poles merge forming branch cuts. One can use the asymptotic expansion 
ip{z) ~ \ogz — l/2z + . . . (valid for |argz| < tv/2) to find the zero-temperature limit of 
( |1.4|) (ignoring the contact terms) : 

G* d ~ log \u 2 - k 2 \ - i ^j^0(^ - k 2 ) sgn u . (1.6) 

For spacelike momenta, \u\ < k, the imaginary part of is exponentially suppressed. 

Even though the retarded Green's function in 4d cannot be found explicitly, the lo- 
cation of its singularities and thus the "dispersion law" of thermal excitations at strong 
coupling can be determined precisely. As shown in |10[ , this amounts to finding the quasi- 



normal frequencies of dilaton's fluctuation in the dual near-extremal black brane back- 
ground as functions of the spatial momentum! . Recently, there has been considerable 
interest (inspired by the AdS/CFT duality conjecture) in studying quasinormal modes of 
Schwarzschild- AdS black holes! [^2],|l^,|lJ,^5|,|l^JT3,|l8l . Here, we consider the non-compact 
case directly. This corresponds to finding quasinormal frequencies for an infinitely large 
AdS black hole. 



Note that the poles (|1.5|) of G^ d (u}, k) coincide with the quasinormal frequencies of BTZ black 



hole [ 11 1 . 



3 References to the early works on quasinormal modes in asymptotically AdS spacetime as well 



as to works considering bulk dimension other than five can be found in [10] , |12|. 
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Computing the scalar quasinormal frequencies in the black brane background is equiv- 
alent to solving the two-parameter singular spectral problem for the equation with four 
regular singularities (Heun equation). In this paper, we solve this problem by using the 
elegant method based on Pincherle's theorem. The eigenvalue equation is written in terms 
of continued fractions and solved numerically. Similar approach had been used by G. Jaffe 
P3|j in 1933 to find the spectrum of the hydrogen molecule ion. Later, it has been applied 



to the problem of quasinormal modes in asymptotically flat spacetimes by E.W. Leaver 

m ■ 

The paper is organized as follows. In the next section, we formulate the problem of 
finding the quasinormal modes of near-extremal black three-brane in terms of the singular 
boundary value problem for Heun equation. In section 3 this boundary value problem 
is solved by analyzing the associated linear difference equation and applying Pincherle's 
theorem. The results (quasinormal frequencies ui n and the "dispersion law" w(k)) are 
presented in section 4. The discussion follows in section 5. 



2. Quasinormal modes and the boundary value problem for Heun equation 

The metric corresponding to the collection of iV parallel non-extremal three-branes in 
the near-horizon limit is given by 




where f(r) = 1 — r 4 /r 4 . According to the gauge theory /gravity correspondence, this 
background with the parameter of non-extremality tq is dual to the M — 4 SU(N) SYM 
at finite temperature T = tq/tiR 2 in the limit N — * oo, gy M N — > oo. 

Using the new coordinate z = 1 — r^/r 2 and the Fourier decomposition 

e -i^-* M z), (2.2) 

the equation for the minimally coupled massless scalar in the background (|2.1| ) reads 

1 + (1-^) 2 , A 2 q 2 

+ z(l -z)(2- z) + 4,2(1 -z)(2- z) 2 * k ~ 4,(1 - z)(2 - z) * h = ° ' (2 ' 3) 

where A = uj/ttT, q = \k\/ixT. Eq. ( |2.3| ) has four regular singularities at z = 0, 1, 2, oo, the 
corresponding pairs of characteristic exponents being respectively {— zA/4, iA/4}; {0,2}; 
{-A/4, A/4}; {0,0}. 
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Quasinormal modes are defined as solutions of Eq. ( |2.3|) obeying the "incoming wave" 
boundary condition at the horizon z = and the vanishing Dirichlet boundary condition 
at spatial infinity z = 1. The first condition singles out the exponent v^p = —i\/4 at 
z = 0. 

The most straightforward way to find quasinormal modes would be to construct a 
local series solution (pi oc (z, A) to Eq.(|2.3|) with the exponent ^q 1 "* near the origin, prove its 
convergence at z = 1, and determine the eigenfrequencies A n (g) by solving the equation 
<£>z oc (l, A) = numerically. This approach works quite well for the low- level eigenfrequen- 
cies, and in fact it has been successfully used in a number of publications on quasinormal 
modes in asymptotically AdS space. Here we would like to solve the above eigenvalue prob- 
lem in a somewhat different way which, in our opinion, is more appealing both analytically 
and numerically. 

By making a transformation of the dependent variable 



(j){z) =Z 4 (z_ 2) 4 y( z ) , 

Eq. ( p73|) can be reduced to the standard form of the Heun equation |2l|.p2 

a(3z — Q 



y" + 



7 5 

- + — 

z z — 



1 + z -2 



y + 



z(z-l)(z-2)- 



0. 



where a, /3, 7, e depend on A, 

A(l + i) 



a 







1 



1 - 



iX 



5 = -l, e=l-A/2. 



(2-4) 



(2.5) 



(2.6) 



(2.7) 

on the interval [0, 1] has 
with A = has received 



4 ' 2 
and Q is the so called "accessory parameter" given in our case by 

„_Q 2 A(l-z) A 2 (2-z) 
V 4 4 8 

We would like to determine values of A and q for which Eq. 
solutions obeying the boundary conditions y(0) = 1, y(l) = 0. 

Before turning to the solution, let us remark that Eq. (|~ 
some attention previously in connection with the so called "glueball mass spectrum" in 
QCD3 [p3[ , [p4],|[25|l, [p6| . Approximate analytic expressions for the "glueball" massess 
squared M% = —q\ were found either via the "brute force" WKB calculation [|26| or by 
appealing to the unpublished results of the RIMS group |2«| (both approaches give sat- 
isfactory agreement with the results obtained by numerical integration). Introducing a 
non-zero A and using the "incoming wave" boundary condition makes the problem consid- 
erably more complicated, as is evident from Eqs. (|2.5|) and 



3. Solving the boundary value problem 

3.1. Local solutions 

The Frobenius set of local solutions near each of the singularities can be easily con- 
structed. 

At z = 0, the local series solution corresponding to the index v = and normalized 
to 1 is given by 

oo 

y (z) = J2 a n(X,k)z n , (3.1) 

n=0 

where oq = 1, di = and the coefficients a n with n > 2 obey the three-term recursion 
relation 



where 



a n+2 + A n (X) a n+1 + B n (X) a n = 0, (3.2) 

= _(n + l)(2^ + e + 3(n + 7 )) + g 
V ! 2(n + 2)(7i + l + 7 ) v ; 

B n (X) = + + , . (3.4) 
nV ; 2(n + 2)(n + l + 7 ) V ; 

The series ( |3.1| ) is absolutely convergent for \z\ < 1 and, in general, is divergent for \z\ > 1. 

The condition for convergence at \z\ = 1 involves parameters of the equation and will be 

investigated below. 

At z = 1, the difference of the exponents has an integer value, and we expect the local 
solution there to contain logarithms. Indeed, the set of local solutions is given by 

yi (z) = (1 - zf (l + 6^(1 -z) + 4 1} (1 - zf + ■ ■ ") , (3-5) 

y 2 (z) = 1 + 6i 2) (l -z) + h yi (z) log (1 - z) + 4 2) (1 - zf + ■ ■ ■ , (3.6) 
where b{ 1] = (q 2 - A 2 + 3A(1 - i))/12, bf ] = (-q 2 + X 2 + A(l - i))/4, h = -(q 2 - A 2 ) 2 /32, 

(1 2) 

and the coefficients b n , n > 2, can be found recursively from the relation similar to 
the one in Q3.2p . Solving the differential equation Eq. ( |2.5| ) essentially means finding the 
connection between the sets of local solutions. 

3.2. The connection problem 

From the general theory of linear differential equations, it follows that three solutions 
yo(z), yi(z), D2{z) are connected on [0, 1] by the linear relation 

yo{z) = A(X, q) yi{z) + B(X, q) y 2 (z) , (3.7) 
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where A, B are independent of z. The vanishing Dirichlet boundary condition for quasi- 
normal modes at z = 1 implies B = 0, which is the equation for eigenfrequencies X n . If 
A = A n , yo(z) is proportional to yi(z), and thus it is simultaneously a Frobenius solution 
about z = (with exponent 0) and a Frobenius solution about z = 1 (with exponent 2). 
Such a solution is called a Heun function |27| . We learn that quasinormal modes of black 



branes are Heun functions. Unfortunately, the connection problem for the Heun equation 
remains unsolved, and explicit expressions for the coefficients A, B are unavailable, with 
the exception of some special cases@. There is, however, an indirect way of determining for 
which values of A and q the connection coefficient B vanishes. This is achieved through an 
analysis of convergence. 

3. 3. The analysis of convergence 



The convergence of the series ( |3.1|) can be analyzed by studying the large n asymptotic 
behavior of the linear difference equation (|3.2|) . One finds! that Eq. ( |3.2|) possesses two 
linearly independent asymptotic solutions of the form 



oo (1) 



n ^ n 

s=0 



4 2) ~ n" 3 



.(2) 



s=0 

(12) 

where the coefficients c s ' can be found recursively using the asymptotic expansion of 
A(n), B(n). In particular, we have c[ 1] = (q 2 - X 2 + A(l + i) + 2A 2 (1 + i))/4, cf } = 
(12 — q 2 + X 2 — 3(1 + i)X)/4. Equation is called a minimal solution to Eq. ( |3T^ ) , and 
Eq. ( |3.9p represents a dominant one. This distinction reflects the property 

a (1) 

lim -4^- = , (3.10) 



4 Schafke and Schmidt [28] give the connection coefficients in terms of the n — > oo limit of a r . 



obeying Eq. ( |3.2j ) . Although the asymptotic behavior of a n near n = oo can be studied by means 
of local analysis, it only determines a n up to a A-dependent coefficient which is essentially the 
very quantity we are looking for. 

5 Necessary information about second order linear difference equations can be found, for ex- 



ample, in the excellent review paper by Wong and Li [29| . 
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and will be useful later on. Using Eqs. (|3.8|) , 



we obtain 



L n+1 



a (1) 



a 



1 

2 

(2) 
n+1 



2 + Re A 
2^ 



+ o 



(2) 
Cln 



1--+0 - 

77, 1 f ' ~ 



(3.11) 



(3.12) 



Equations (|3llD , (|37l2D imply! that generically the series (|3.1|) converges absolutely for 
|z| < 1 for any value of A. This proves that the equation 



J2°<n(\,k) =0 



(3.13) 



71 = 



can indeed be used to compute the quasinormal frequencies numerically. What is more 
interesting, however, is that in some cases the radius of convergence increases. It happens 
when the minimal solution (|3.8| ) exists. In that case the series ( |3.1| ) converges absolutely 
for \z\ < 2 and thus represents an analytic function (the Heun function) in that region. 

Thus, even though the connection coefficient £>(A, q) in Eq. ( |3.7|) remains unknown, 
finding zeros of £>(A, q) is equivalent to finding a condition under which the minimal solution 
to Eq. ( |3.2|) exists. Such a condition is conveniently supplied by Pincherle's theorem |U 
which states that the minimal solution exists if and only if the continued fraction 

So (A) Si (A) S 2 (A) 



A (X)- A 1 (X)- A 2 (X)- 
converges. Moreover, in case of convergence one has 

o-n+i S n (A) S n+1 (A) S n+2 (A) 



(3.14) 



a, 



A n (X)— A n+ i(A)- A n+2 (X)- 

The right hand side of Eq. ( |3.15|) is generated by a simple algorithm. Define r r 
Then Eq.( |3.2| ) can be written as 

S n (A) 

T — — 

A n (X) + r n+1 ' 

Applying ( |3.16[ ) repeatedly, one gets Eq. ( |3.15| ) . Now, setting n = in Eq. ( |3.15| ) we have 

Q S (A) Si (A) S 2 (A) 



(3.15) 

dn+l/ 0-n- 

(3.16) 



2-iX A {X)- A X {X)- A 2 (X)- 



(3.17) 



See go | , §2.37. 
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This is the transcendental eigenvalue equation which determines the quasinormal frequen- 
cies. Equation ( |3.17| ) can be solved numerically with great efficiency using a variety of 
methods (see |32| , which also includes a discussion of the error analysis). Here we use a 
nonlinear backward recursion which amounts to breaking the continued fraction ( |3.16|) by 
setting r n to zero@ for some large n* and computing backward to get tq. Stability can be 
checked by choosing a larger value of and repeating the calculation. 



4. Quasinormal frequencies 

For q = 0, the lowest 15 quasinormal frequencies A n obtained by solving Eq. (|3.17| ) 
numerically are listed in Table 1 and shown in Fig. 1. The results suggest that the number 
of frequencies is infinite, and that their large n asymptotic behavior is given by the simple 
formula 

A± = A±±2n(l T z), (4.1) 
where A^ ~ ±1.2139 — 0.7775 i. In terms of the original variable u, Eq. ( |4.1| ) is written as 

w± =a;±±27rTn(lT*), (4.2) 
where the coefficient in front of the parentheses is the bosonic Matsubara frequency and 

The lowest eigenfrequency in Table 1 can be compared with the result of Horowitz 
and Hubeny [12] for a large 5d Schwarzschild- AdS black hole. Normalizing the r_|_ = 100 



entry in Table 1 of p| appropriately (A = u;R 2 /r + , R = 1), we have A^ 100 ) « 3.119627 



2.746655 % which is fairly close to our result A (oo) « 3.119452 - 2.746667 i. 

The dependence of the lowest five quasinormal frequencies on q for q G [0, 100] is 
shown in Figs. 2-4. All five branches stay above the line Re A = q (Re uj = |k|) slowly 
approaching it in the low-temperature limit q — > oo. The imaginary part of the branches 
tends to zero in the same limit. Figures 2-4 bear resemblance to the dispersion law of 
thermal excitations in a weakly coupled Yang-Mills plasmaR 



7 To improve the convergence of continued fractions, one can instead set r n to its asymptotic 
value at n*, r n = 1/2 — (2 + A)/4n* + ■ • •, as suggested by Nollert |33| . 

8 Note, however, that in this paper we consider correlators of the gauge-invariant operators 
rather than the gluon propagator. 
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n 


Re X n 


Im A n 


1 


±3 


119452 


-2 


746676 


2 


±5 


169521 


-4 


763570 


3 


±7 


187931 


-6 


769565 

1 \J *J \J KJ \J 


4 


±9 


197199 


-8 


772481 


5 


±11 


202676 


—10 


774162 


6 


±13 


206247 


-12 


775239 

1 l KJ U %J iJ 


7 


±15 


208736 


-14 


775979 


8 


±17 


210558 


—16 


776515 


9 


±19 


211943 


—18 


776919 


10 


±21 


213025 


-20 


777232 


11 


±23 


213896 


-22 


777489 


12 


±25 


213896 


-24 


777489 


13 


±27 


213896 


-26 


777489 


14 


±29 


213896 


-28 


777489 


15 


±31 


213896 


-30 


777489 



Table 1: The lowest quasinormal frequencies A n for q = 0. 
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Fig. 1: The lowest 15 quasinormal frequencies in the complex A-plane for 



All quasinormal frequencies have negative imaginary parts. This fact can be proven 
rigorously using an argument based on the "energy-type" integral [|34J or by writing 
Eq. (|2.3[) in the Regge- Wheeler form and essentially repeating the proofs for AdS black 



holes given recently in [12], [14|. It reflects the stability of the near-extremal metric against 



9 



5 10 15 20 25 

Fig. 2: Re A vs q (with the interval Aq = 0.5) for the lowest five quasinormal 
frequencies. The zeros approach the line Re A = q as q — ► oo. 




a scalar perturbation. On the field theory side, the negative sign corresponds to the damp- 
ing of the plasma excitations. 

The computed frequencies appear to be symmetric with respect to the imaginary A 
axis, i.e., the eigenvalue equation ( |3.17| ) seems to give complex-conjugate pairs of solutions 



in terms of the variable iX. This symmetry does not seem to be explicit in Eq. ( |3.17| ) 



9 The symmetry, however, is obvious in the original form of the equation ( |2.3|) . I thank 
A. V. Shchepetilov for pointing this out to me. 
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5. Discussion of results 

The results are compatible with expectations outlined in the Introduction. The re- 
tarded Green's function has infinitely many poles in the complex w-plane whose location 
depends on the spatial momentum. There is a "mass gap": for small enough values of A 
the propagator G R (iv, k) is analytic in agreement with Eq. (|1.3|) . In the low-temperature 
limit A —>■ oo, q — > oo singularities merge, forming branch cuts of the zero-temperature 
propagator ( |1.2| ) as in the BTZ case. 

The distribution of quasinormal frequencies in the complex ui plane for an asymptoti- 
cally AdS background appears to follow a much simpler pattern than the one corresponding 
to the asymptotically flat case (compare Fig. 1 in [[20| ). In particular, Chandrasekhar's 
"algebraically special" solution is absent: there are no frequencies with Re A = 0. Obvi- 
ously, all statements about the behaviour of higher order modes are conjectural. It would 
be very desirable to confirm the asymptotic formula (|]T) analytically, possibly by using 



the complex WKB method. We remark, however, that the analogous problem remains 
unsolved even in the much studied case of a Schwarzschild black hole in flat space. 

Spectra of excitations in more realistic theories can be similarly studied provided 
their gravitational duals are known explicitly. It would also be interesting to determine 
the poles of the current and energy-momentum tensor correlators at finite temperature, 



10 The zero-temperature limit ( |1.2| ) of G R (u, k) can be obtained from the Heun equation using 
the Langer-Olver asymptotic expansion (see [|To| ). 
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whose hydro dynamic limit has recently been computed in |J . 

Finally, studying gravitational quasinormal modes may be important in investigating 
the stability of the non extremal black brane background. 
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